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I review some recent results on hadronic decays of charmed mesons. These decays illuminate a wide range of physics, including the 
absolute normalizations of b and c decays, the understanding of D° mixing, the coupling of heavy mesons to hadronic decay products, 
and light hadron spectroscopy. 

1. INTRODUCTION 

The lightest open charm mesons — the D°, D + , and D + s — decay through the weak interaction, and the majority of their 
decays are to final states containing only hadrons. These mesons and their hadronic decays lie at an intersection of weak 
and strong physics: the easiest detection of charmed states for flavor physics is through hadronic decays with large rates and 
simple topologies, while the decay processes themselves provide important information on hadronic spectroscopy and strong 
interactions. 

Four groups of measurements are discussed here. Two of them, absolute branching fractions of D mesons and strong phases 
in hadronic decays, are important for allowing precision studies of the CKM matrix and other physics to reduce systematic 
uncertainties. The other two, the studies of two- and three-body D decay amplitudes, shed light on how long distance hadronic 
physics affects the visible results of short-distance weak processes. 



2. EXPERIMENTS 

Recent activity in the field has been dominated by the charm experiment CLEO-c and the B-factories BaBar and Belle. 
Detailed understanding of hadronic charm decays places a premium on efficient and pure particle identification, high resolution 
on daughter hadron momenta, well-understood initial states, and large datasets, which these experiments can provide. 

The CLEO-c experiment at the CESR symmetric e + e~ storage ring took data in the center of mass energy range between 3.67 
and 4.26 GeV to investigate topics in both open charm and charmonium physics. Primary data collection in the charm region 
occurred between 2003 and 2008. Of particular relevance here are 818 pb _1 of </f(3770) data for D° and D + studies at threshold 
and around 600 pb _1 in the region of 4.17 GeV collected for D + s physics. This yielded approximately 3 million D () D , 2.4 
million D + D 1 , and six hundred thousand D^D* pairs. By running just above threshold for the production of various species, 
it could be guaranteed that D°, D + , or D + s mesons would always be produced with their antiparticles. This is one of the most 
prominent features of the CLEO-c dataset and enables powerful partial reconstruction (or "tagging") techniques, as the typical 
efficiency for full reconstruction of a single D meson (in one of a menu of clean decay modes) is between 5 and 10% and the 
presence of the antiparticle is implied for those events. In addition, the low final state multiplicity of D and continuum light 
quark production processes at this energy means combinatoric background is generally small. 

The BaBar and Belle experiments are detectors at the PEP-II and KEKB asymmetric e + e~ facilities, respectively. Both 
experiments began taking data in 1999; Belle is ongoing, while BaBar's run concluded in 2008. Although the datasets of the 
experiments are collected in the BB threshold region, there is an appreciable rate for charm production, both from B decays 
and from continuum e + e~ — > cc events. The analyses mentioned here use datasets ranging from 210 to 550 fb _1 . The large 
luminosities compensate for the relatively lower efficiencies compared to CLEO-c, and statistical precision can be comparable 
or better at the B-factories. The boost of the center of mass frame at these experiments also allows them to reduce backgrounds 
by taking advantage of displaced vertices due to D meson lifetimes. 
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Figure 1: Measurements of the branching fractions of three reference D decay modes: (top left) D° — > K~n + , (top right) D + — > K~n + n + , 
(bottom) D + s -> K-K + n + . Ref. (D is "CLEO 07"; Ref. (2) is "CLEO 08"; Ref. is "BaBar 08"; and Ref. (4) is "Belle 07." The comparison 
is made to the PDG 2004 average [5 1 as that is the last edition not to include new input from CLEO, BaBar, or Belle. 

Table I: Branching fractions for three important reference decays of D mesons, showing the evolution of between 2004 and 2008. All 
branching fractions are in percent. 





S{D° -» K~n + ) S(D + K-7t + tt + ) 


S(D+ -» K K + k + ) 


PDG 2004 fit |H 


3.80 ±0.09 9.2 ±0.6 


4.4 ± 1.2 


CLEO-c 


3.891 ± 0.035 ± 0.069 1 1 1 9.14 ± 0.10 ± 0.17 □ 


5.50 ±0.23 ±0.16 (2) 


BaBar 


4.007 ± 0.037 ± 0.072 (3) — 




Belle 




4.0 ± 0.4 ± 0.4 H (prel.) 


PDG 2008 fit 


3.89 ±0.05 9.22 ±0.21 


5.50 ±0.28 



3. ABSOLUTE BRANCHING FRACTIONS 

Certain high rate, low multiplicity decays of D mesons serve as standard calibrations for decays of charm and bottom quarks. 
With very few exceptions, all measurements of branching fractions or decay widths for D or B mesons depend eventually on 
these reference decays. For example, measurements of \V c b\ via exclusive B — > D^tv decays require knowing the absolute 
width for this process, which requires precision knowledge of D° and D + decay rates. 

Over the last few years there has been significant improvement in our knowledge of these rates, as illustrated in Fig. [T] and 
Table [I] Between the 2004 and 2008 editions of the Review of Particle Physics El 151. the uncertainties on the three reference 
branching fractions shown in Table [I] have been improved by factors of 2 to 4. 

CLEO-c obtains yields of "single tag" (a D or D decay is reconstructed in a given mode) and "double tag" (both a D and a 
D are found in specified modes) events, and compares single and double tag yields for various decays to extract the branching 
fractions. For N decay modes, yields are obtained for 2N single tag decays (D and D are considered separately) and for A^ 2 
double tag decays. The yields can be predicted in terms of branching fractions, the total number of DD pairs, and the 
detection efficiencies, the last being determined from Monte Carlo simulations. There are A^ 2 + 2Af measurements and N + \ 
parameters, providing an overconstrained and well-conditioned set of equations which are solved for the branching fractions 
using either^ 2 minimization or likelihood maximization. This general technique has been used for D° and D + at 3.77 GeV/c 2 
|Q] and for D + s at 4. 17 GeV/c 2 [2| with three, five, and eight decay modes respectively. 

BaBar measures T>(D° — > K~n + ) using a partial reconstruction technique |3| with the high-rate process B° — > D* + (X)f~v (X 
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represents possible extra particles that are not reconstructed). The soft pion from the D* + decay is reconstructed, as well as the 
lepton f . The D* + momentum is estimated using the fact that the energy release in D* + — > D°n + is very small so there is a 
strong correlation between the momenta of the pion and the D* + . The inclusive yield of events is determined from a fit to the 
implied neutrino mass, assuming the initial B was at rest. Signal D° — > K~n + candidates are then reconstructed in these events. 
Backgrounds are obtained using below-resonance data for continuum contributions and Monte Carlo simulation of BB decays; 
the background expectations are verified using events where the soft pion and the lepton candidates have the same sign instead 
of opposite charges. The result is S(D° -> K-jt + ) = (4.007 + 0.037 ± 0.072)%. 

Belle has a preliminary result for S(Z)+ — > K~K + n + ) that uses a double-partial -reconstruction technique (4). They exploit 

continuum production of D* + D~ 1 (2536) (along with the charge conjugate), with D* + — > D + S y and D~. —* D K. Four-momentum 

— * 

conservation can be used to infer the presence of the D s or D if the other three particles are reconstructed. The observed 
yield of D Xy(Dj)inferred and jD^K(D inferred depend on the branching fractions of D and D+, respectively, and since the 
total yield of D* + D~ 1 (2536) is the same in both cases the ratio gives information on Z)+ branching fractions. The result is 
S(D+ -> K-K + n + ) = (4.0 ± 0.4 ± 0.4)%. 

As is evident from TablejlJ at the moment measurements of D° and D + branching fractions are systematics-limited (generally 
by determination of tracking and particle identification efficiencies) while measurements are still statistically limited. To 
round off the set of important charm hadron reference branching fractions, it would be very desirable to see a result using similar 
techniques for S(A t t — > pK~n + ), for which no model-independent measurements currently exist. 

3.1. Normalization of D' Decays 

The older BaBar results on &{D + S — > (fm + ) ||7][8] have so far not been mentioned. These measurements impose specific mass 
requirements on the kaon pair in the K~K + n + final state. Because there is a non-negligible contribution of other processes to the 
K K + n + signal in this mass region (notably — > /o(980)7r + ) the choice of mass window will affect the obtained result. These 
results are hard to relate to those using the full K~K + n + phase space because, in particular, the ratio of branching fractions will 
not be the fit fraction of D + s — > <pn + obtained from an amplitude analysis. It is worth noting that the 2008 edition of the PDG no 
longer couples S(D^ — > <fm + ) to &(D+ — > K~K + n + ) via a fit fraction for this reason. 

For many experimental purposes, however, it is desirable to restrict the allowed phase space to the <p region to improve the 
ratio of signal to background. There are two proposals to relate branching fractions in restricted and full phase space regions. 

One approach was taken in Ref. HI, which gave, in addition to the full branching fraction, "partial" branching fractions 
which are obtained by requiring events to have K K + mass within fixed windows around m$. Except as necessary to derive 
reconstruction efficiencies, no assumptions are made about the resonant structure in this region — in particular the measured 
branching fraction is not corrected for S((p — > K K + ) and the scalar contribution is fully included. 

The other, longer term approach is to perform careful amplitude analyses to produce a precise description of the shape of the 
D + s — > K~K + n + distribution over the entire Dalitz plot. Other experiments can then implement this description in physics event 
generators such as EvtGen j9), and determine the efficiency of their specific selection choices relative to the full phase space 
by Monte Carlo simulation. Towards this end BaBar [ 10] and CLEO are both studying the — > K~K + n + Dalitz plot and 
have reported preliminary results. Projections of their fits are shown in Fig. [2] and the fit fractions listed in Table |H| The results 
agree in general terms, although there is some tension in the region of the (f> due to the scalar contributions (in particular BaBar 
quotes a large systematic uncertainty in the fit fraction due to /o(980)7r + ). For this purpose the critical question is how much 
the two descriptions of the Dalitz plot differ on the rate at different points in phase space, which can only be checked when full 
amplitude and phase information are made available. 

4. STRONG PHASES FROM QUANTUM CORRELATIONS 

One of the primary techniques for exploring D° mixing is to look at the time evolution of the "wrong-sign" decay D°(t) — > 
K + n~, where D°(t) indicates a meson, which was produced in the definite flavor eigenstate D°, after proper time t has elapsed. 
This amplitude receives two contributions: 
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Figure 2: Projections of data and amplitude fits onto Dalitz varables for the process D* — > K K + n + from preliminary CLEO (top) and BaBar 
(bottom) 1 10] analyses. The CLEO analysis uses 14 x 10 3 events, and the BaBar result uses 67 x 10 3 decays. 



Table II: Fit fractions of various decay modes in preliminary D* — > K~K + tc + Dalitz analyses from BaBar [ 10] and CLEO-c. All fit fractions 
are in percent. The uncertainties on the BaBar results are statistical and systematic, respectively; those on the CLEO-c results are statistical 
only. 

Mode BaBar CLEO-c 

K*(S92)°K + 48.7 ± 0.2 ± 1.6 47.4 ± 1.5 
£*(1430)° K + 2.0 ± 0.2 ± 3.3 3.9 ± 0.5 
£•(1430)° K + 0.17 ±0.05 ±0.3 — 
/o(980)7r+ 35 + 1 + 14 28.2 ± 1.9 
cfm + 37.9 ± 1.8 42.2 ±1.6 

/o(1370)tt + 6.3 ±0.6 ±4.8 4.3 ± 0.6 
/ (1710>r + 2.0 ±0.1 ±1.0 3.4 ±0.5 
/ 2 (1270K 0.18 ±0.03 ±0.4 — 



• the direct but doubly-Cabibbo-suppressed (DCS) decay D° — > K + n~ , 

• mixing followed by the Cabibbo-favored (CF) decay D° — > K + n~. 

In the limit of small mixing the first amplitude remains near its initial value, while the second grows linearly; this differing 
dependence permits a measurement of mixing parameters from the shape of the D°(t) — > K + n~ rate. Because of the additional 
(and larger) amplitude from the DCS decay, the mixing term is amplified by interference and becomes measurable at first order 
in the mixing parameters. However the size of the interference depends on the relative phase of the two amplitudes as well as 
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their magnitude, and relating this measurement to the fundamental mixing parameters requires knowing this phase. The weak 
contribution to the phase is known, leaving the QCD component to be determined. 

This can be done by exploiting the quantum correlations between D pairs at CLEO-c, as detailed in Ref. IfTTI . The process 
e + e~ — > i^(3770) produces a state with the vector quantum numbers J PC = 1 . When this decays to the pseudoscalar pair 
D°D , the daughters are produced in a L — 1 angular momentum state, which implies that the product of the daughter CP states 
must be odd. We can rotate the basis and think of the ip (3110) decay as producing a pair of mass eigenstate mesons D\D2, 
which in the absence of CP violation are also CP eigenstates: 

|P°) ± [P°) 
V2 

One can apply the standard tagging techniques in this new basis. If a CP eigenstate decay is tagged, this projects out the 
opposite CP component on the other side — which is a sum or difference of flavor eigenstates, whose decay amplitudes will 
interfere. In particular the resulting total amplitudes (and rates) to common D° and D final states, for example K + n ~, will be 
dependent on the phase between the two contributing amplitudes. 

The CLEO-c analysis |[T2l to extract the strong phase 5 between D° — > K + n~ and D — > K + n~ involves reconstructions 
of three classes of D° decays: CP eigenstates (e.g. K'K* for CP = +, KQtt for CP = -), semileptonic decays (which are 

unambiguous flavor tags as the lepton has the same charge sign as the parent charm quark), and K^n*. The CP correlation 

n —° 

modifies the time-integrated double tag yields from what is expected in the "uncorrelated" case where the D and D decay 
independently. For example, CP = + versus CP — + is suppressed to zero, while D° -> K 7r + vs. D — > Xe v e is modified 
by a factor (1 - ry cos 6 - rx sin S), where 6 is as above, r is the magnitude of the amplitude ratio between DCS and CF K + n~ 
decays, and x and y are mixing parameters representing the mass and width difference of the mass eigenstates. 

In principle enough information is present in this analysis to determine all the input parameters (r, 5, and the mixing pa- 
rameters x and y, as well as relevant branching fractions). In practice CLEO-c event yields are insufficient to do this, so 
external inputs are used. In the standard fit the branching fractions as well as Rm = {x 2 + y 2 )/2 and the time-integrated ratio 
R ws = Y(D° -> K + tt-)IY(D -> K + n-) are fixed to precise determinations from the B-factories and CLEO-c. In addition, in 
the standard fit x sin 6 is fixed to zero as it is extremely weakly determined and sin 6 is expected to be near zero (the associated 
systematic uncertainty on cos 6 is +0.03). 

The likelihood contours from the standard fit are shown in Figure [3] The likelihood function is peaked at cos 6 ~ 1, but has 
a significant fraction of the area outside the physical region. In addition, the uncertainties are not Gaussian. By projecting the 
likelihood onto the cos 5 axis and integrating 95% of the area in the physical region, CLEO-c obtains cos 6 > 0.07 or |<5| < 75° 
at 95% confidence level. This analysis provides the first measurement of 6. 

An "extended" fit can also be performed incorporating external information on r and the mixing parameters x and y. This 
significantly reduces the uncertainty on y and permits them to float xsind. At 95% confidence level CLEO-c finds xsind € 
[0.002, 0.014] and 5 e [-7°, +61°]. This result depends on external measurements that can be expected to become more precise, 
and hence is less useful than the standard fit for incorporation into averages. 



5. DECAY MECHANISM STUDIES 

Leptonic and semileptonic decays of heavy quarks are much better understood theoretically than hadronic decays, because 
the relevant weak interactions are demonstrably short distance processes and the daughter leptons do not participate in the strong 
interaction. The effects of the strong interactions can thus be absorbed in parameters such as the decay constant for leptonic 
decays and form factors for semileptonic decays, which can be predicted using various techniques. 

In comparison, hadronic decays involve at least four quarks. These can arrange themselves into hadrons in many ways; in 
particular the produced hadrons can rescatter into other states. These "long distance" processes occur at the hadronic scale 
instead of the electroweak scale, and are much less theoretically tractable than short distance interactions and subsequent two- 
quark behavior. 



5 



34 th International Conference on High Energy Physics, Philadelphia, 2008 



3950108-008 




Figure 3: Likelihood function for determination of the K ~n + strong phase via quantum correlations at CLEO-c 1121 . Plot a) shows the 
likelihood projected onto the cos 5 axis; the hatched region shows 95% of the area in the physical region (| cos 6\ < 1). Plot b) shows the two 
dimensional contours in y and cos 6; the lines are spaced lcr apart, where cr = -\/Ax 2 - The result cos <5 « 1 is stable under excursions in y, but 
the uncertainties are not Gaussian. 



Long distance effects can be large. For example, the short distance contribution to D mixing is highly suppressed by the GIM 
mechanism, but long distance effects due to flavor SU(3) violation can induce significant mixing and are likely the dominant 
cause of the observed lifetime difference of the two mass eigenstates fOl . 



5.1 . Observation of £>+ -> pn 

The process D + s — > ph is the only kinematically allowed decay of a ground state charmed meson to baryons and provides a 
unique probe of hadron dynamics in weak annihilation decays Ifl4l [I5l [T6l . The recently observed large branching fraction of 
order 10~ 3 , in contradiction to initial predictions of order 10~ 6 , raises the question of whether final state rescattering dominates 
short distance annihilation, or whether applications of the PCAC hypothesis ignoring states heavier than the pion fail in this 
case. 

This decay has been seen for the first time at CLEO-c ifTTl . No attempt is made to find the neutron or antineutron from the D + s 
decay; instead its presence is inferred from kinematic variables, in particular the missing mass opposite the rest of the D^D* 
event, which is fully reconstructed. This is analogous to the procedure used for the CLEO-c measurement of T>{D + S — > p + v fl ) 
|[T8l and many of the methods are shared. A tag D s is reconstructed in a fully hadronic mode, and combined with a photon 
candidate to see if the missing mass of that combination corresponds to mass of the as far unreconstructed D + s . If so, and in 
addition the photon is kinematically compatible with arising from the decay of the D* + or D*~, a proton track is searched for in 
the drift chamber; particle identification is done using specific ionization. If found, the missing mass opposite the D~yp system 
is evaluated. Kinematic fits to the assumed intermediate particles are used to improve resolution and reject background. The 
results are shown in Fig.|4j there are thirteen events in the signal region of 900-980 MeV/c 2 . Sidebands in the tag D s mass are 
used to evaluate backgrounds; no sideband events fall in the signal missing mass region. Thirteen events and no background 
yields a branching fraction T>{D + S — > ph) = 1.30 + 0.36^q [q x 10~ 3 , which is appreciable, especially given the phase space 
suppression. 



5.2. Two-Body Decay Amplitude Analysis 

Although the inclusion of long distance effects means that short distance Feynman diagrams for weak decays of charmed 
quarks should not be taken seriously for computational purposes, they are still useful for visualizing the flow of quark flavor 
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Figure 4: Observation of D* — » ph at CLEO-c 1 17]. Left: missing mass squared of tag D s y combination, showing peak at m 2 D for true 
£>**DJ events. Candidates within the red arrows are then investigated for D + s — > ph candidates. Right: missing mass opposite D s yp system 
showing peak at m„. The solid line shows the expected signal distribution. 




Figure 5: The six flavor topology diagrams used to analyze D — > PP data in Ref. [22]: (top left) tree, (top center) color-suppressed tree, (top 
right) annihilation, (bottom left) annihilation with singlet emission, (bottom center) exchange, (bottom right) exchange with singlet emission. 

during a decay in the "flavor topology" approach. Coupled with flavor SU(3) symmetry, this can be a powerful tool for 
understanding the relationships between different hadronic decay rates. 

CLEO-c has recently assembled a complete picture of Cabibbo-favored and singly-Cabibbo-suppressed open charm decays to 
pairs of pseudoscalar mesons H] El OS El ED; many of these are first observations, especially for D + s decays. These branching 
fractions have been analyzed by Bhattacharya and Rosner [22] in a comprehensive framework assuming SU(3) symmetry to 
relate D°, D + , and D + s decays. Six flavor topologies were assumed, and are detailed in Fig. [5] The results are shown in Fig. [6] 
Among the conclusions are: 

• the T diagram has the largest magnitude, but C and E are also fairly large: \C\ ~ 0.75|7"|, \E\ ~ 0.6|r|; 

• there are large phases between T, C, and E (e.g. 150° between C and T); 

• the annihilation diagram is significantly smaller and approximately antiparallel to E: A 0.3E; 

• there are two allowed solutions for SE, both of which are small compared to E; 

• there is one solution for SA, which is of comparable magnitude to A. 
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Figure 6: Results of analysis of D — > PP decays by Bhattacharya and Rosner 1221 . Left: complex amplitude for T, C, E, and A amplitudes 
(see text and Fig.|5]for explanation of amplitudes). Center: construction to obtain S E amplitude; allowed solutions lie at the intersection of 
all circles. Two solutions are evident. Right: similar plot for SA; there is only one solution. 



It should be emphasized that these conclusions are only valid for two-pseudoscalar decays, and that because of long distance 
effects these should not be interpreted as statements about short distance physics. (For example, the large E amplitude does not 
imply that short distance W exchange amplitude is also large.) 



6. DALITZ PLOT STUDIES 

The decays of D mesons provide interesting information on light hadron properties. All possible routes to the same final 
state interfere, so phase information is available. Because the initial state is a pseudoscalar, the angular distributions are highly 
constrained. The case of decays to three pseudoscalars is particularly popular to study ("Dalitz analysis"), as there are only 
two independent kinematic variables and no information hidden in the final state daughter spins. In general it is assumed that 
the variation of the amplitude for the decay X — > ABC can be decomposed into a the coherent sum of amplitudes for a few 
intermediate states, where two of the daughters arise from a resonance and the third is a spectator (for example, — > n + n + ji~ 
could have contributions from D + s — > /o(980)7r + , /2(1270)7r + , and p(770)7r + ); this is termed the isobar model. In particular 
all amplitude and phase variation is assumed to come from the intermediate resonance which is generally parametrized using 
Breit-Wigner or Flatte functions. 

This description, while very useful, has recently been recognized as inadequate, especially for describing the S -wave com- 
ponent of decay amplitudes. In particular a coherent sum of interfering Breit-Wigner or Flatte lineshapes does not conserve 
unitarity, and the resonances in the scalar sector are so wide that this is a major issue. There is also a debate on the nature of 
the scalar sector below 1 GeV/c 2 , where structures referred to as the cr and k are suggested to contribute to im and Kn scatter- 
ing, respectively. Their importance in charm decays is unclear. In the absence of more sophisticated techniques, the S -waves 
in D decays are typically parametrized as constant ("non-resonant") amplitudes with ad hoc "resonances" added where the 
non-resonant description fails drastically. 

One technique for unitarizing the S -wave amplitudes, called the Zf-matrix approach, has been successfully applied to charm 
decays in the past. It relies on outside information on phases and amplitudes of low-energy scattering to derive amplitudes, 
and is explicitly unitary. Unfortunately this method is complex, the physical meaning of the results is usually not especially 
transparent, and it requires external input — the experimental result is not self-contained. 

Recently another approach, reliant on high statistics, has been pioneered by the E791 Collaboration Il23l : the "model- 
independent partial wave analysis" (MIPWA). If two amplitudes with different angular dependence interfere, and a model 
(Breit-Wigner, etc.) is available for the phase and magnitude of one, the complex amplitude for the other can determined, up to 
a twofold phase ambiguity, without assumptions on its energy dependence. (Requiring smooth evolution of the amplitude with 
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Figure 7: Scalar K~n + (left) and 7r + 7r + (right) amplitudes in the process D + — > K~n + n + , as obtained by CLEO-c 1241 . For the A" - 7r + plot, 
the points indicate the MIPWA results for the amplitude excluding the ^(1430) contribution; the dot-dashed line shows the A*q(1430) Breit- 
Wigner shape; the solid line indicates the sum of the MIPWA and A^(1430) amplitudes; the dotted line shows the result of an isobar model 
fit including a k and a constant nonresonant K^n* amplitude; the dashed line shows the k component of that fit. For the n + n + plot, the points 
show the result of the MIPWA fit for the 7r + 7r + amplitude; the dashed line shows the analytic amplitude used in the MIPWA fit for the K~n + 
amplitude; the dotted lines show ±lcr variations on the analytic amplitude used in the isobar model fit to the K~n + amplitude. 



energy can help choose between solutions.) Vector and tensor resonances in D decays are generally modeled well by Breit- 
Wigner lineshapes with well-known masses and widths, which can be fixed in order to determine the amplitude and phase of 
the S -wave components. Ambiguities can still arise in such an approach but the fact that the raw amplitudes are made available 
makes this a very powerful technique. 

Here I mention two such analyses. The first is a CLEO-c study of the D + — > K ~n + n + decay ll24l using 141 x 10 3 events. The 
dominant contributor is found to be the K~n + S -wave. The standard MIPWA is modified by explicitly including a K (143Q) 7i + 
contribution with a Breit-Wigner shape for the ^f (1430)°; this is done so that the remaining scalar amplitude varies slowly. 
This amplitude is observed to be large at low K~n + mass; however it does not follow the expectation for a k alone, and can only 
be modeled in that framework by adding a large non-resonant component that interferes significantly with the k at low mass. 
From this it is unclear that the k is a compelling model for the K ~n + S -wave in this decay. 

The fit also requires an appreciable n + n + S -wave component (fit fraction 10-15%), whose parametrization is taken from fits 
to pion scattering data and which has significant variation with n + 7T + invariant mass. The MIPWA technique was used to obtain 
amplitudes for this component as a cross-check of the nominal analytic function, and the two show good qualitative agreement. 
The results for the K ~n + and n + n + S -wave amplitudes are shown in Fig. [7] 

The second is a preliminary BaBar analysis of the — > n^n + ii~ decay using 13 x 10 3 events, discussed in more detail in 
Ref. Il25l . This decay is dominated (83% fit fraction) by the tt + t:~ 5-wave. For the these studies are particularly compelling 
because of the large fraction of hadronic decays where the Cabibbo-favored ss pair is not manifest, where one mechanism to 
explain this would involve long distance ss — > uu, dd through states where those mix, such as the scalar sector. The only large 
non-scalar contribution is from /2(1270)7r + . 

The results of the MIPWA are shown in Fig. [8] There is a clear peak for the /o(980) with the expected variation in phase for 
such a dominant resonance. In the 1400 MeV/c 2 region, there is also significant activity. However, the low mass amplitude 
or phase variations are not drastic, and careful study is needed to gauge the compatibility of these results with significant <x 
production. 
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Figure 8: Scalar ;r + 7r amplitude in the decay — > 7r + ;r + 7r from a preliminary BaBar analysis [25 1, showing large / (980) contribution and 
structure in the 1.4 GeV/c 2 region. 



7. SUMMARY AND OUTLOOK 

In the last few years, the large datasets that have become available at CLEO-c, BaBar, and Belle have enabled sophisticated 
studies of hadronic D decays to be performed. These results have reduced branching fraction uncertainties by factors of 2 to 
4, provided first measurements of strong interaction phases between two body decays, improved our understanding of decay 
dynamics, and shed light on light hadrons. There is more to learn in all these areas. The high-luminosity BES-III experiment 
has recently started taking data and will have much to contribute to our understanding of the hadronic decays of open charm. 
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